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Abstract. The dense cores which conceive and cradle young stars can be explored through continuum emission from associated 
dust grains. We have performed a wide field survey for dust sources at 1.2 millimetres in the p Ophiuchi molecular cloud, 
covering more than 1 square degree in an unbiased fashion. We detect a number of previously unknown sources, ranging from 
extended cores over compact, starless cores to envelopes surrounding young stellar objects of Class 0, Class I, and Class II type. 
We analyse the mass distribution, spatial distribution and the potential equilibrium of the cores. For the inner regions, the survey 
results are consistent with the findings of previous narrower surveys. The core mass function resembles the stellar initial mass 
function, with the core mass function shifted by a factor of two to higher masses (for the chosen opacity and temperature). In 
addition, we find no statistical variation in the core mass function between the crowded inner regions and those in more isolated 
fields except for the absence of the most massive cores in the extended cloud. The inner region contains compacter cores. This 
is interpreted as due to a medium of higher mean pressure although strong pressure variations are evident in each region. The 
cores display a hierarchical spatial distribution with no preferred separation scale length. However, the frequency distribution 
of nearest neighbours displays two peaks, one of which at 5000 AU can be the result of core fragmentation. The orientations 
of the major axes of cores are consistent with an isotropic distribution. In contrast, the relative orientations of core pairs are 
preferentially in the NW-SE direction on all separation scales. These results are consistent with core production and evolution 
in a turbulent environment. Finally, we report the discovery of a new, low-mass Class object candidate and its CO outflow. 
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1. Introduction 

As one of the nearest star-forming regions, the p Ophiuchi 
molecular cloud complex has been the target of numer- 
ous investigations. A prime focus has been the densest part 
of the L 1688 cloud which harbours a large number of 
young stellar objects, studi ed at near- t o mid-infrared wave- 
lengths JVrba et all Il975t lEliasI Il978l IWilking et alJ Il989t 
Bars onv e t al.| |l989t JComeron et al. 1993; Luhman & Rieke 
1999 ; IWilking et alJEoOll iBontemps et alJl200lHAllen et alj 
2002) and at (s ub)millim etre wavelengths (Andre et al. 
1993| lAndre & Montmerld 1 19941 iNiirnberger et alJ 1 19981 
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IWilson et al]ll999fe Tachiha raetal]l2000l) . To investigate how 
these stars are conceived, we wish to to relate their properties to 
those of the embedding cloud of molecular gas and dust which 
both nurture and obscure the stars. However, technical capa- 
bilities restrict the field over which area-covering surveys can 
be undertaken, bearing the danger of picking peculiar objects 
prevalent in the "more interesting" parts, thus biasing the re- 
sults obtained from statistical investigations. In order to further 
exclude such selection effects, we have performed an unbiased 
survey of a wide area around the L 1688 molecular cloud for 
dust continuum sources, including the dense molecular cores as 
well as the areas which apparently do not harbour dense molec- 
ular material. 

The entire star formation complex extends over several de- 
grees on the sky (Loren 1989a), containing a few major clouds. 
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The L 1688 cloud is situated about one degree south of the 
p Ophiuchus star itself. This 'main cloud' of p Oph covers 
an area of roughly 480 arcmin 2 a nd has been dissecte d into 
about a dozen cloud components jMezgeret alJll992l) . This 
area was surveyed at 1.3 mm bv lMotte"etalT i Il998l hereafter 
MAN98), who uncovered 62 starless cores and 41 circumstellar 
structures. These observations had a resolution of just 1 1 arc- 
sec, corresponding to 1,400 AU. The inferred distribution of 
masses of these cores was found to be comparable to the stel- 
lar initial mass function, suggesting that stellar masses are de- 
termined at conception. Subsequently, the results of a larger, 
som ewhat more sensitive survey at 0.85 mm were presented 
bv ljohnstone et al] l l2000l) (hereafter J00). This survey covered 
~ 700 arcmin 2 with a resolution of 14 arcsec, identifying 55 
cores. This survey was recently followed up by a much more 
extensi ve (4 square degrees) b ut significantly more shallow sur- 
vey bv ljohnstone et al] J2004I) (hereafter J04), mainly extend- 
ing to the north/north-west of the main cloud. 

In the present 1.2 mm survey, we cover an area of 
4,600 arcmin 2 . This constitutes a significant extension to the 
above works although the resolution is ~24". Our immediate 
goals are to (i) derive the core masses (©, (ii) estimate the 
mass relative to a critical Bonnor-Ebert sphere or Jeans mass 
by including the core sizes (©, (iii) quantify the spatial distri- 
bution and orientations (|5] & fBJ, and (iv) quantify the mass 
distribution ({j6j. In addition, we report the discovery of a new 
low-mass candidat e Class source and its CO outflow (^}. A 
companion paper ( Khanzadv an et all20 04) presents the results 
of a near-infrared imaging search for protostellar H2 outflows 
over a large section of the present survey area. Finally, the core, 
outflow and protostellar properties will be related in detail in a 
forthcoming paper. 

Quoted distances to the p Ophiuchi cloud lie in the r ange 
D _ 1 25 - 165 pc, as rece ntly summarised by iRebull et al] 
j2004 . lde~Geus et all JT989) analysed photometry of the stars 
in the nearby Sco OB2 association. It was shown that the 
Ophiuchus dark clouds are on the near side of Upper Scorpius, 
at ~ 125 pc. Hipparcos data confirm the distance to Upper 
Scorpius and yield a distance of ~ 128 pc to the quadruple p 
Oph system Jde Zeeuw et alJll999l) . Here, we shall follow the 
suggestion by Reb ull et al] (|2004) and take the weighted aver- 
age distance of 130pc. 

In the following, we will designate all well-defined low 
mass structures as 'cores' (with mass below about 50 M Q ) and 
apply the term 'clump' to more massive structures defined 
through low resolution CO mapping. A clump usually corre- 
sponds spatially to a cluster of cores. 

2. Observations 

The observations were carried out using the SIMBA bolometer 
array at the SEST telescope on La Silla/Chile during an observ- 
ing run lasting from 2002 July 7 to 12. SIMBA is a 37-channel 
bolometer array, operating at 1 .2 mm, and provides a HPB W 
of ~24" (as measured on maps of Uranus). The weather was 
good and fairly stable with zenith opacities ranging between 
0.17 and ~0.3. In total, 78 maps in the fast mapping mode 
were taken (i.e., no wobbler was used). Typical map sizes were 



1200"x800" with a scanning velocity of 80"/sec and 8" steps 
between subscans. A number of maps were taken with a scan- 
ning velocity of 160"/sec (size 1800"x800") in order to im- 
prove sensitivity for extended structures. Pointing and focus 
were checked regularly; pointing corrections were typically a 
few arcsec and always below 10". Sky dips were taken every 1 
to 2 hours, and Uranus was mapped several times for calibra- 
tion purposes. 

The data reduction was performed with MOPSI, a soft- 
ware package developed by Robert Zylka, following standard 
bolometer data reduction principles. First, noisy or dead chan- 
nels were identified and deleted. Then, a low order baseline was 
subtracted from the data and an initial de-spiking was made. 
The data direct from the telescope consist of the true signal con- 
volved with th e response of the electronics, hence requirin g de- 
convolution jReichertz et aljEoOlt IWeferling et al]l2002i) . We 
chose to reconstruct the signal at a specific time using only data 
from a given window (typically 6 seconds) around this time, in- 
cluding the data recorded during the "turning" of the telescope 
between two on-map sub-scans (in the standard deconvolution 
only data taken during the on-map subscans are used, and the 
signal is reconstructed from the entire on-map sub-scan). This 
deconvolution was found to provide significantly better results 
than the standard deconvolution, for two reasons: the restric- 
tion to a shorter time span reduces noise that goes into the re- 
construction of a given data point, and including the time dur- 
ing the turning of the telescope ensures a proper reconstruction 
of the data at the edges of the map, as sky fluctuations during 
turning (which are convolved into the data recorded during the 
on-map subscans) are properly taken out of the data recorded 
during the on-map subscans. Then, after the deconvolution, the 
data corresponding to off-map data taken during the turning 
of the telescope were deleted. After further low-order baseline 
corrections, correction for atmospheric extinction and correc- 
tion for gain-elevation effects, correlated sky brightness vari- 
ations (sky-noise) were removed. Here, an iterative procedure 
was used: for each iteration, we used the resulting mosaic of the 
previous iteration as an input model for the brightness distribu- 
tion. Following another careful de-spiking and baseline sub- 
traction, the data were then combined into a mosaic, weighting 
individual maps according to their noise (1/rms 2 ). After each it- 
eration, for each sub-map the polygons defining emission-free 
regions to be used for baseline fitting were carefully checked 
and, if necessary, adjusted, followed by a check of the degree 
of the polynomial used for baseline subtraction. 

The noise level in the final mosaic is of the order of 
10 mJy/beam throughout most of the map. Close to the edges, 
it degrades to 12-15 mJy/beam. For a point source, assuming 
a flux-to-mass conversion as outlined below, this corresponds 
to a 3<t detection limit of the order of 0.014 M G . In terms 
of column density, the lcr noise level corresponds to about 
1.3xl0 21 cirT 2 . This is slightly better than the 1.3 mm survey 
done by MAN98, but somewhat less sensitive than the J00 
850yum survey; it is significantly more sensitive than the recent 
J04 850/vm survey. 

Similar to other (sub)millimetre surveys, our new map is 
sensitive only to sufficiently small structures (see also J00 for 
a discussion of this issue). In the case of SIMBA fast mapping 
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Fig. 1. SIMBA 1.2 mm continuum mosaic image of the p-Ophiuchi cloud. Intensity scaling is in mJy/beam (linear). Prominent 
clumps are marked as A, Bl/1, C, D, E, and F. Dashed lines mark some filamentary structures east of the major clumps, with three 
more compact sources (MMS038, MMS051, and MMS126) within them also labeled. MMS126 is identified as a new, low-mass 
Class object candidate in Sect.0 



observations, this is due to two effects: the subtraction of the 
sky emission removes any flux due to large scale, more or less 
uniform emission, and a high-pass filter applied to the signal 
from the bolometer will further suppress large-scale emission 
features. On top of these technical effects, baseline subtraction 
may further contribute; however, baseline subtraction was done 
very carefully in order to minimise this. Filtering of low fre- 
quencies and residuals from the deconvolution of the data are 
responsible for residual negative map values close to the bright- 
est parts of the mosaic (dipping to — 20 mJy/beam (-2 <x) in the 
most negative patches), a feature generally found in SIMBA 
maps. 



3. Data Analysis 

The final 1.2 mm mosaic is shown in Fig.^ In li ne with pre- 
vious (sub)millimetre maps, a large variety of sources is seen, 
ranging from extended cores over more compact features to un- 
resolved point sources. In those regions which were also cov- 
ered by the MAN98 and J00 surveys we generally see a very 
good agreement between the different data sets (limited by the 
poorer angular resolution of our new data). 

In the part of the mosaic which had not been previously ob- 
served, besides a number of fainter, extended structures, two 
more remarkable features turned up. On the one hand, there 
is a filamentary structure (dashed lines in Fig. ^ to the east 
of the previously revealed cores Oph-B/C/E/F, containing two 
more compact sources (MMS126 and MMS038/MMS051). 
MMS126 appears to be a newly discovered, low-mass Class 
source (see Sect.Qfor details). 

The second remarkable structure is the presence of a bright 
dust core MMS041 at the tip of the cometary dark cloud 
L 1696B. A similar core was already found in the neighbour- 
ing cloud L 1696 by MAN98. 

3.1. Source identification 

The variety in source morphologies makes a detailed source 
identification mandatory: e.g., using simple source detection 
algorithms with a certain detection threshold will not iden- 
tify extended structures, which are clearly present, but might 
have maximum surface brightnesses of the order of the l<x 
noise level. To account for this, sources were identified using a 
wavelet decomposition similar to the method used by MAN98. 



The final map was split into 5 planes using the 
TRANSFORM/WAVE and EXTRACT/WAVE tasks provided 
by the ESO/MIDAS wavelet context package. We will refer to 
the most compact scale as scale 1, and scale 2, 3, 4, and 5 to 
the next more extended scales. 

We then used the 2-D adaption of the clumpfind algorithm 
developed by Will iams et alJ (1 1994b to identify features in the 
5 planes individually. The lower limits for source identifica- 
tions were set to 3-4 times the rms fluctuation of emission-free 
parts of the respective planes, with contours spaced by about 
1 rms. The minimum number of pixels required to identify a 
sufficiently bright feature as a source was set according to the 
characteristic scale of each plane (i.e., to very small for scale 
1, and getting bigger for larger scales). The thus detected fea- 
tures were then visually inspected to reject false detections due 
to residual scanning artifacts (found to be restricted to scales 1 
and 2) or other misidentifications. 

From the list of features identified at the various scales a 
final source list was compiled, including all features identified 
in scale 1, and all sources in the more extended scales which did 
not overlap with sources identified in the next more compact 
scale (i.e., for each feature we determined the part of the clump 
area as given by the clumpfind routine, in which the intensity 
was greater than 1/2 the maximum intensity, and checked if it 
overlapped with the corresponding area of all features of the 
next more compact scale). This resulted in a list of 139 sources 
(MMS001 toMMS139). 

To this list we added 4 sources, MMS140 to MMS143, af- 
ter visual inspection: while MMS141 is just too faint to be 
detectable in the individual planes, MMS140, MMS142, and 
MMS143 are small cores adjacent to large bright features, 
which make them hard to detect as taking out the extended 
scales creates negative haloes around large bright cores. 

A full finding chart showing the location of all 143 sources 
is shown in Fig. IA.1I (online version only), and closeups of 
crowded regions are shown in Figs. IA.2llA.3l and lA. 41 (online 
version only). 

Table |A~T1 (online version only) lists the 1.2 mm sources 
found in our survey in the order they were identified by the 
wavelet decomposition plus clumpfind technique. In addition 
to the source position, it holds the integrated source flux and 
the peak flux determined in various ways, the major and mi- 
nor axis and the position angle of the major axis, cross- 
identifications to the sourc es identified by Mo tte et al.l ([1998) 
and Johnstone et al. ( 2000), and a comment on association with 
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Fig. 2. Model of the 1.2 mm continuum emission derived from the superposition of models for the 143 sources identified from 
the SIMBA data. Intensity scaling is in mJy/beam (linear). 
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Fig. 3. Residual after subtracting the model source distribution (Fig.|2j from the 1.2 mm mosaic (Fig.[0. Intensity scaling is in 
mJy/beam (linear). 



known YSOs. In Section |A] we give notes on the individual 
sources. 

3.2. Measuring source properties 
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Fig. 4. The locations of all the 143 cores detected here (squares) 
compared to all the cores detected by J00 (asterisks). The coor- 
dinates are centred on -16:26:58.4, -24:45:36 (2000). A, B, and 
E mark the location of prominent clumps. 

We then performed a detailed modelling of the overall 
brightness distribution using Gaussian sources, the result of 
which is shown in Fig. [2] (Fig. [5] shows the residual, after sub- 
tracting the model source distribution from the mosaic). In 
some cases, more than one Gaussian component centered on 
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Fig. 5. The locations of the 111 starless cores detected here 
(squares) compared to the the starless cores detected by 
MAN98 (asterisks). The circle is that employed to separate in- 
ner cores from the outer regions. 



the source position had to be assumed in order to properly re- 
produce the shape of the sources (e.g., a bright source with a 
radial intensity distribution falling off as a power law will not 
be satisfactorily fitted by a single Gaussian). 

The sources we detect within the inner core-crowded region 
are generally coincident with those found in previous surveys. 
Figs.|4]and|5]show that this is true for separate samples of all 
cores and starless cores. In particular, we recover all the cores 
from J00, despite the different observing wavelength. This con- 
firms the reliability of our core-finding technique. 
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The fluxes of the sources were determined by first subtract- 
ing the model for all other sources from the mosaic, and then 
integrating over a certain area around the source within some 
given radius or (for more complex sources) within a certain 
polygon around the centre of the source. Similarly, sizes and 
position angles for each source were determined by (single) 
Gaussian fits to the mosaic with the model for all other sources 
subtracted. 

Peak fluxes were calculated in four differing ways and are 
also given in Table lA.il the column labelled as F pea k (a) gives 
the peak flux obtained from searching the maximum pixel value 
in an ellipse (with the minor and major axes as given in the 
table) around the peak position; Fp ea k (b) is the same as (a), 
but using only the sum of scales 5 down to the scale where the 
source was detected first; F pea k (c) gives the maximum pixel 
value in the same ellipse after re-binning according to the size 
of the sources minor axis (i.e., by FWHM m j n /24"pixels); F pea k 
(d) gives the maximum pixel value of the model of the source. 
As done in deriving the total source flux, the models of all other 
sources were subtracted from the mosaic before determining 
the respective peak fluxes. 

Depending on the source morphology and brightness, one 
or the other method can be expected to yield more accurate 
results. Relatively bright, compact sources should be best rep- 
resented by Fp ea k (a), whereas faint, extended sources should be 
best measured by F pea k (b) or (c). Generally, the results agree 
well, apart from an obvious overestimate for large, low surface- 
brightness sources using the first method and systematic un- 
derestimates for compact sources using the second and third 
methods. 

At this point it might be worth adding a note on the de- 
pendency of the effective mass/column density sensitivity limit 
on source size. The surveys of J00 and J04 identified sources 
setting a fixed threshold in peak flux, i.e., column density. In 
contrast, we also include sources whose surface brightness is 
nominally lower than, e.g., 3cr, but which are still clearly recog- 
nised because they are extended (the MAN98 analysis using a 
wavelet decomposition might be comparable in this respect). 
This introduces a size-dependent column density limit: for a 
source which has twice the radius R, we can rebin the mosaic 
by 2x2 pixels, effectively reducing the noise by a factor of two. 
I.e., the column density detection limit scales as R _1 , which 
implies a mass detection limit scaling as R 1 (whereas a fixed 
surface brightness cutoff implies a mass detection limit scaling 
as R 2 ; always assuming that all sources have the same temper- 
ature and dust properties). In this sense, we can expect that our 
method of source identification yields, for a given mass, a more 
complete census than applying simply a column density cutoff, 
because it also includes the more extended, hence fainter (in 
surface brightness), sources. 

Our method of source identification yielded significantly 
more sources than identified by J00 in areas covered by both 
surveys, despite the somewhat poorer angular resolution and 
slightly lower sensitivity. This is due to two reasons. First, J00 
explicitly filter out more extended emission before applying 
clumpfind, thus removing extended features and lowering the 
signal also for more compact features. Second, they limit their 
identification to features with a given column density threshold, 



again biasing against low surface brightness, extended sources. 
Detailed comparison of our data with the J00 map before filter- 
ing out the extended emission shows that virtually all features 
that we identify are also seen on the J00 map, even if not iden- 
tified in their list of objects. We have marked those features as 
'n.i.' in the 'J00' column of Tab. [A~T1 

4. Core Properties 

4.1. Masses and column densities 

The mass of each core is estimated on assuming it to be pro- 
portional to the measured 1 .2 mm flux. The radiation is taken to 
be emitted from cool dust grains and the intervening medium 
is optically thin. The standard formulae for black bodies then 
yield the core masses from the fluxes Si. 2mm expressed in 
Janskys (e.g MAN98 and J00): 



M r = 0.58 x S 



1 :2mm 



K\ .2mm 



0.005 cm 2 ; 
/12K 

exp Nr 



1 



D 



130 pc 



Mo 



(1) 



Here, we take a fixed opacity at 1.2 mm of *• = 0.005 cm 2 g -1 
and a fixed dust temperature of Tj = 20 K. It should be noted 
that these quantities may vary across the cloud by a factor of 
order two, as fully discussed by MAN98 and J00. 

The derived masses are plotted as a function of core size 
in Fig. [6] Here, we do not attempt to deconvolve the appar- 
ent cloud dimensions and define the mean observed clump an- 
gular radius of 6 = 0.5 (FWHM„ MJ x FWHM m ,>,) 1/2 from the 
FWHMs of the major and minor axes given in Table lA. fl 

Although the inferred masses and radii are widely dis- 
tributed, a mass-radius correlation is evident in Fig. [6] 
However, it should be remarked that large, low-surface bright- 
ness features might have been supressed by the observational 
procedure and data reduction, as discussed above. Hence, the 
apparent correlation is possibly only produced by a deficiency 
of extended, moderately massive features. (It should be noted 
that the smallest core radii presented in the diagram are at 
the limit of our resolution and subject to significant error). At 
small radii and masses, the lower envelope to the distribution 
of sources reasonably well reproduces the size dependency of 
the mass detection limit as argued above, Mi; m oc R. To in- 
vestigate a possible mass-radius correlation, locii of constant 
average column density, corresponding to logM oc 21ogR e , 
are overplotted in Fig. [6] This demonstrates that at low radii 
and masses the detection limit yields a lower envelope to the 
distribution which influences the correlation. 

Splitting the sample into cores within the inner, crowded 
region (marked by crosses, corresponding to the region inves- 
tigated by MAN98 and J00) and those spread over the ex- 
tended cloud (diamonds) hints at distinct lower envelopes for 
cores in the inner and outer regions, corresponding to aver- 
age atomic hydrogen columns of 8 x 10 21 cirT 2 (inner region) 
and 4 x 10 21 cirT 2 (outer regions). We believe that this is at 
least partially due to observational limitations, as low surface- 
brighntess features are hard to identify in the crowded central 
region. 
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Fig. 6. Upper Panel: The mass and radius derived from the 
observations for each core. The crosses correspond to the 79 
cores found within the central circle of radius 0.2° shown in 
Fig. El an d diamonds to those 64 outside. Lines of constant 
column density are overplotted for reference as dashed lines, 
with columns of atomic hydrogen from 3 x 10 21 crrT 2 (right) 
to 3 x 10 22 crrT 2 (left). The solid line marks the detection limit. 
Lower Panel: The masses and radii for just the starless cores, 
consisting of 62 within the circle (crosses) and 49 outside (dia- 
monds). 



The number distribution for the core columns is displayed 
in Figure^ This demonstrates that the offset in column density 
between the inner and outer regions is real and extends to the 
entire distribution and not just the lower bound to the column 
density. It is not caused only by a (possibly obervationally in- 
troduced) deficiency of low column density cores in the inner 
region, but also by the absence of high column density cores in 
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Fig. 7. The number distribution of the mean column density of 
the cores (also proportional to the surface pressure predicted 
if a sub-sample were critical Bonnor-Ebert spheres, as given 
from Eq. |3j. Upper Panel: The full sample of 143 cores, with 
79 within the inner crowded core and 64 external to it; lower 
Panel: the same distribution for just the starless cores (111 in 
total, 62 within the inner crowded core and 49 external to it). 



the outer region, which would have been easily detected. The 
inner cores possess a systematically higher column density by 
a factor of approximately two on average. 

This suggests that a form of segregation, present in stel- 
lar clusters, is also taking place in the clumps in p Ophiuchus. 
However, as we will demonstrate in Section [6] the mass dis- 
tribution is very similar in the inner and outer regions. The 
segregation is not in the mass but in the column density or 
pressure, hence compactness of the cores. That is, the cores 
must be somewhat more compact to survive in the crowded re- 
gions. This could also suggest that cores, in general, cannot be 
treated as isolated objects although the compacter cores, which 
exist sufficiently long, may have reached a relatively indepen- 
dent state. This scenario , however, must be reconciled with the 
dynamical time scales JBelloche et all200ll) which allow little 
time for core-core interactions. 

The existence of two exceptionally massive, high density 
cores in the inner region can be considered to be the result of 
physical processes rather than core overlap along the line of 
sight. As can be seen from Fig. |6] these two cores are both 
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massive and distinctly compact and so lie in a distinct region of 
the mass-radius diagram. These are the cores MMS001 (known 
as SMI and SM1N) and MMS004 (SM2) within the most 
crowded Clump A. T hey appear to out line a cavity around the 
bright, young star SI (Grasda len et al]l973l) . This suggests that 
they consist of swept-up and compressed gas, locally triggered 
into the collapse process. 

Note also that restricting the analysis to just the starless 
cores (Figs. [6] and lower panels) does not alter the above 
results. However, the majority of the very low mass compact 
cores, both within and outside the crowded region, already con- 
tain young stellar objects. This suggests that such low-pressure 
cores only form or survive by being gravitationally bound to a 
point source. 

4.2. Core Support 

Are the cores ephemeral or in a prolonged state, close to equi- 
librium? Some cores could be internally supported by a uni- 
form thermal, turbulent or magnetic pressure and contained by 
a roughly equal external pressure. In such two-phase media, 
the ambient medium may be quite diffuse, warm atomic gas. 
However, this cannot apply to the majority of the cores since 
the internal density generally rises significantly towards their 
centres (MAN98). Nevertheless, an ambient medium of con- 
stant pressure could be maintaining a few of the diffuser low- 
mass cores in the outer regions. With a limited variation in core 
temperatures, uniform pressure would roughly lead to a rela- 
tionship of the form logM oc 31ogR c , i.e., steeper than that 
corresponding to cores of constant column. This is clearly not 
observed. 

Alternatively, the derived masses can be compared to those 
expected from isothermal spheres of the same size which 
are in hydrostatic equilibrium. In particular, there is a maxi- 
mum cloud mass which is stable to perturbations in pressure 

■ II I I II I 

(Ebert 1955; Bonnor 1956). The maximum stable mass of a 
Bonnor-Ebert sphere, which assumes self-gravitating, isother- 
mal cores bounded by an external surface pressure, is Mbe 
1.18 cr 4 /(G^t 2 P l J 2 ) where P s is the surface pressure and cr is 
the one-dimensional vel ocity dispersion, h ere equivalent to the 
isothermal sound speed dGalli et al.l2002l) . Evaluating in terms 
of a typical molecular cloud pressure, we obtain 



M»e = 1.2 



20 K 



2 10 -10 dyne cm 



-1/2 



Mo. (2) 



In terms of size, the equilibrium sphere will be stable if the 
radius exceeds the critical radius Rbe = 0.41 G Mbe/o 2 - 
Assuming a fixed core temperature (MAN), we thus predict a 
shallower relationship log M oc log R f if the cores are close to 
the critical state. 

In a similar fashion, we could assume that the surface pres- 
sure is quite uniform and allow the temperature of the cores to 
vary by a small factor, consistent with Eq.|2] Then, eliminating 
cr we obtain the critical radius 

1-1/4 



Rbe = 5.1 x 10 J 



M BE ' 


1/2 




[1 M Q _ 




2 10 



Ps 



1-10 



dyne cm 



AU, 



which predicts the relationship log M oc 2 log R c . 

The data, displayed in Fig. [6] indicate a wide range in core 
masses for any given core size. Under the strict assumption that 
all cores are critical Bonner-Ebert spheres, have the same tem- 
perature and only thermal internal motions, this would imply a 
range in surface pressures exceeding 100. The loci of constant 
mean column presented in Fig. [6] also correspond to surface 
pressures of 2 x 10~ 10 dyne cirT 2 (upper left line), 2 x 10" 
dyne cirT 2 (middle line) and 2 x 10" 12 dyne cirT 2 (lower right 
line). In contrast, J00 found a much narrower mass-radius cor- 
relation for cores in the inner region, consistent with a smaller 
range in surface pressure. 

The core mass is significantly less sensitive to the core size 
than that predicted by a constant column. Least squares fits of 
the form 



log(M/M ) = a log(R c /10, 000 AU) + b 



(4) 



(3) 



yield a = 1.11 ± 0.18 and b = -0.16 ± 0.06 for all starless 
cores (correlation coefficient is 0.50) with a — 1.45 ± 0.22 and 
b - -0.24 + 0.05 for the outer starless cores (correlation coeffi- 
cient is 0.70) and a = 1.49 ±0.30 and Z> = 0.07 + 0.12for the in- 
ner starless cores (correlation coefficient is 0.51). Figure[S]dis- 
plays the mass-radius relation for just the outer starless cores. 
Obviously, the correlation is not very stringent, as also indi- 
cated by the fairly low correlation coefficients. The offset be- 
tween the power laws for the inner and outer cores corresponds 
to a difference in size by a factor of ~ 1.2-1. 3, which using 
Equn. translates to a pressure difference of a factor of ~2- 
3. This is significantly less than the surface pressure variations 
which can be inferred from the range of core masses at some 
given radius or as have been derived by J00. Random pressure 
variations within the inner and outer cloud gas seem to domi- 
nate the systematic increase in pressure when going from the 
outer to the inner part, as one could expect for a medium being 
subject to random, turbulent motions. 

Although the observational limit tends to flatten the mass- 
radius relation, we cannot exclude an underlying intrinsic law 
of the form logM oc logR f , as given by critical Bonnor-Ebert 
spheres with a constant velocity dispersion. In this case, how- 
ever, the wide range in masses implies a range of 100 in the sur- 
face pressures to maintain the cores in equilibrium, as given by 
equation |3 This variation in pressure could be present, result- 
ing from highly supersonic turbulent motions (see §[SJ. In fact 
more likely would be a steeper relation than indicated by the fit, 
approaching log M oc 2 log R e , as expected for critical Bonnor- 
Ebert spheres immersed in a constant external pressure medium 
having varying velocity dispersions; due to the cr 4 dependence 
of Mbe only moderate variations would be needed to create the 
more than 1 order of magnitude dispersion of mass at a given 
radius. Again, this would more likely be a sign for turbulent in 
addition to thermal motions, as it appears unlikely that starless 
cores exhibit large enough temperature variations. Moreover, 
in a turbulent medium, cores will be continually forming and 
dispersing, and may temporarily resemble bound cores. In re- 
cent gravo-turbulent simulations, starless cores were generally 
gravitationally unbound, suggesting that gravitat ional collapse 
occu rs promptly after gravity becomes dominant JKlessen et alJ 
120051) . 
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Fig. 8. The mass-radius relation for the starless cores in the 
outer regions. The solid line is the least squares fit and the 
dashed line is suspected to be the detection limit with M <x R c . 



5. Spatial Distribution 

The identified cores are strongly clustered at first sight. To 
quantify this, we employ a two-point correlation function treat- 
ing each core as a point object in space. Following J00, we 
determine the number of core pairs, i/ ( /(r), with separation be- 
tween log(r) and log(r) + d\og{r). This is compared to the pre- 
dicted distribution for a random sample of cores spread over 
the apparent volume, H,(r). The two-point correlation function 
is then defined as 

H d (r) 



H r (r) 



L. 



(5) 



For the random sample, we consider a random surface dis- 
tribution across a sheet of the survey size and shape. Fig.0 
shows that the correlation function is positive on scales under 
~ 10 5 AU. For comparison, the data obtained by J00 are also 
plotted on the Figure. 

The survey size is seen to influence the correlation with the 
present survey displaying more large scale structure while the 
J00 survey detected more smaller cores, as expected from the 
observational constraints. J00 found that a power law of the 
form <i> cc Ar" 75 fitted the data well out to 3 x 10 4 AU. The 
new survey suggests that this relationship can be extended to 
above 10 5 AU. Therefore, in contrast to J00, we exclude the ex- 
istence of a preferential scale near to 3 x 10 4 AU, which could 
correspond to a Jeans length. Instead, the most obvious possi- 
ble cause of hierarchical clustering is turbulence. Supersonic 
turbulence would also generate a wide range in pressure and 
column density, as required to interpret the mass-size distribu- 
tion (see §[S}. 

The inferred power law index of -0.63 is somewhat flatter 
than found by J00, although within the error bars except on the 
largest comparable scales. Note that, as remarked by J00, this 
value is close to the value measured for galactic clustering of 
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Fig. 9. The two point correlation function as defined in the text 
for the p Ophiuchus cores. The error bars correspond to VN 
statistics. The dotted straight line corresponds to a power-law 
of the form a> cc Ar~ 063 . 



-0.668 for scales under 1 ° dMaddox et al.ll99ol) . Furthermore, 
the function O is closely related to the mean surface density 
of companions (MSDC), being approximately proportional for 
values of <£ exceeding unity (Simonl ll997l) . Thus, it is inter- 
esting that the power law radial dependence of the MSDC for 

pre-main-sequence stars in Taurus-Aurigae is 0.6 for scales 

r > 0.04 pc (and corresponds to a fractal point distribution with 
index 1.4). The lower limit separates the bo und systems from 
the un bound stars, and was suggested by iLuhman & Riekd 
(1999) to correspond to both a Jeans length and the size of 
bound molecular cores. 

For samples of Ophiuchi stars, ISimonl (1997) found an 
MSDC index of -0.5 + 0.2 above a break at 5000 AU. On 
the other hand. iNakaiima etaD dl998l) found -0.36 ± 0.06. As 
noted bv lBate etalT ( [l998l) . flattening of the correlation is ex- 
pected on a timescale of 10 5 yr, as unbound stars separate and 
the system expands. 

The frequency of separation between neighbouring cores 
complements the two-point correlation function since it con- 
tains information on all orders of the correlation functions. 
The mean neighbour separation might also be in terprete d as a 
typica l Jeans fragmentation length jGomez et alJl993llLarsonl 
1995). A length scale could be apparent which might be hid- 
den in the distribution involving all pairs. In fact, the mean 
separation between neighbours in our starless core sample is 
17,400 AU and Fig. 1101 indicates that there is indeed a fre- 
quency peak near this separation distance. For comparison, a 
random sample would have a mean separation of 22,100 AU 
(we simulated 500 sets of randomly distributed starless cores 
spread through the entire region). The random sample also in- 
dicates that we cannot interpret the peak in the separation fre- 
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quency as being derived from som physical process which gen- 
erates a break (knee) in a power-law functions with a steep de- 
cline at large distances since a turnover is expected even in ran- 
dom data. 

Of much more significance is the second frequency peak at 
low separations (5000 AU) prominent in both panels of Fig. 1101 
appearing as if the hierarchical core clustering is modified by 
the fragmentation of some close cores. In fact, this separation 
distance corresponds well to the size of the cores themselves 
(see Fig. [(J. Hence, it is not surprising that spatial structure is 
present on this scale. 




3.5 4.0 4.5 

log (Separation in AU) 



5.0 




3.5 4.0 4.5 

log (Separation in AU) 



5.0 



Fig. 10. The number distribution of separations between neigh- 
bouring cores for all cores (upper panel) and the starless cores 
(lower panel). The error bars correspond to -\JN statistics. The 
superimposed lines correspond to the predicted mean num- 
bers for a random sample. The mean separation values are 
15,400 AU (all cores) and 17,200 AU (starless cores). 



6. Mass Distribution 

The distribution of core masses is displayed in Fig.^2f° r a ll 
the starless cores in the sample. We also display in Fig. El the 
mass functions for the inner (62) and outer (49) starless sam- 
ples, defined by a circle of radius 0.2°, located so as to encom- 
pass the crowded inner region studied by MAN98 and J00. We 
find that a similar mass function to that found by MAN98 and 
J00 for the inner zone applies to the entire region. 
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Fig. 11. The number distribution by mass of the 111 starless 
cores in p Ophiuchus. The inset displays their locations. The 
dotted lines are the power law fits of the form dN/dM oc M a . 
Breaks in the power laws occur around 0. 1-0.3 M Q and around 
0.5-1 M . The Salpeter mass function would have a slope of 
a = -2.3 in this format. The error bars correspond to V^V 
counting statistics. The inset displays the core positions in a 
wide region stretching 1° x 1.2°, with the locations of the 
clumps A, B and E marked. 

The core mass function appears to be hard to fit with a sin- 
gle power law, being steeper at higher masses and flattening 
continously towards lower masses. Power law fits with one or 
two breaks are reqired for a satisfying fit. The two-break power 
law fit as displayed in Figs.^2 ancl Elwith two breaks is com- 
parable to the appr oximation for the stellar initial mass function 
(IMF) derived hv lKrounallkOOll) . where the power law indices 
are -0.3, -1.3, and -2.3, although with considerable potential 
variation. Note that the flatter power law holds for core masses 
below ~ 0. 14 M , in comparison to the value of 0.08 M Q quoted 
for the stellar IMF (lKroupal200ll) . 

If we suppose that the break at low core masses is real and 
indeed occurs at around 0.14 M , then the shift to a 0.08 M 
break for the stellar mass function has to be explained. On the 
one hand, it could be due to errors in the assumptions which go 
into the core mass determination. Although a closer distance 
would lower the derived core masses, we have already assumed 
a quite close distance to p Ophiuchus. A larger opacity would 
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Fig. 12. The number distribution by mass and location of the 
cores in p Ophiuchus for 111 starless cores. The dotted circle 
drawn in the inset divides the cores into an inner zone of 62 
objects (diamond symbols) and an outer area of 49 objects (tri- 
angles). The dotted lines, error bars and inset are as described 
in the Fig. II llcaption. 



also reduce the core masses: an opacity of 0.009 cm 2 g -1 would 
move the first breakpoint to 0.08 M for the distance of 130pc. 
However, a somewhat lower dust temperature than the assumed 
20 K in the fiux-to-mass conversion would tend to shift the 
break point to higher masses. Alternatively, if we assume the 
first break and the opacity are accurate, we can speculate that 
just about one half of a core mass ends up constructing the star. 
The other half is dispersed in (i) jets (up to 30%), (ii) dispersal 
by jet impact, (iii) early stellar and disc winds and (iv) other 
stars, brown dwarves and planets. 

Bontemps et al.ld200ll) infer a 2 component power law IMF 
for Class II young stellar objects in p Oph from ISOCAM data. 
The break the IMF is seen at around 0.55 M and separates 
power law slopes with indices -1.35±0.25 and -2.7 below and 
above the break, respectively. This IMF is 'statistically indistin- 
guishable' from the core mass function derived by MAN98 and 
is also consistent with the present data sets. 

Broken power law fits to the core mass function are not 
very well constrained, concerning the number of breaks as well 
as their location. A more continously changing distribution 

such as a loa-normal function might do a better job. In fact, 

[ n 1 j 1 ° J 
Ballesteros-Paredes et al. (2005) derive core mass distributions 

from numerical models of turbulent molecular clouds, finding 
mass functions similar to log-normals. As their technique of 
clump extraction is similar to our approach (taking all density 
enhancements, whether bound or not) both studies should be 
well comparable. Finding a core mass function similar to the 
ones derived from simulations further supports the idea of tur- 
bulence as the main agent in shaping the cloud. 



Whatever the exact, detailed shape of the core mass func- 
tion is, its general, overall behaviour resembles the stellar mass 
function. The MAN98 and J00 surveys both concentrate on 
cores which are argued to be likely gravitationally bound. Our 
survey includes more extended, lower surface brightness fea- 
tures, which are less likely to be bound. However, numer- 
ical simulations showed that the structure of transient den- 
sity peaks forming and dispersing in a turbulent cloud may 
closely resemble hydrostatic, g ravitationally bound objects 
(Ball esteros-Paredes et alJEo 03). so even the cores identified 
by J00 and MAN98 might turn out not to be bound. Any 
scenario which closely correlates the core and stellar masses 
makes presumptions concerning the stratification of the en- 
velopes and their subsequent evolution (i.e. whether certain 
cores might collapse while others disperse). To select cores ac- 
cording to such criteria requires a consistent means of interpret- 
ing the data from cores of very contrasting sizes and masses. 
Moreover, as pointed out by J00, such surveys are inevitably in- 
complete despite our unbiased strategy. Given all these caveats, 
the resemblance of core and stellar mass functions is even more 
remarkable. 

Possibly the most important result of our study is that no 
significant difference is measured in the core mass function be- 
tween the inner and the outer regions where the functions over- 
lap. As shown in Fig. [21 the few highest mass cores lie in the 
inner region. Thus, while the compactness of a core varies in 
space, the mass distribution does not. This must be reconciled 
with the fact that the inferred surface pressures are, on average, 
lower in the external region which would imply a higher Jeans 
mass (neglecting possible temperature differences). 



7. MMS126: a low-mass Class candidate and its 
CO outflow 

MMS126 is a compact but resolved source forming part of a 
roughly east-west oriented, filamentary structure. It appears to 
coincide with a faint, cold IRAS source (IRAS 16253-2429: 
F12: <0.25 Jy; F25: <0.42 Jy; F60: 2.91 Jy; F100: 16.34 Jy). 
Close inspection of the HIRES processed IRAS maps reveals a 
faint source at 25yum, which however coincides with the nearby 
MMS060 = WS B 60. This latter source is also detected by ISO 
at 7 and 14 pm dBontemps et aklEoOll) but no clear source is 
seen at the position of MMS126. Finally, M MS126 appears to 
be dr iving a collimated near-IR H2 outflow dKhanzadvan et all 
2004) and a molecular (CO) outflow (see Sect. l7.1l below). 

Taking these arguments together (apparently cold FIR to 
millimetre source, no clear counterpart at wavelengths shorter 
than 25pm, and association with an outflow) we conclude that 
MMS126 is a very young stellar object, possibly still in the 
Class phase. MMS126 would then be the second Class 
object found in the p -Oph main cloud besides the prototypi- 
cal source VLA1623 dAndre et al.ll993() but is of significantly 
smaller mass: using Eqn.^ we derive a circumstell ar dust mass 
of abo ut 0.2 M Q (VLA1623: 0.6 M , this work: lAndre et alJ 
dl 9931) 1. 
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Table 1. Physical properties of the MMS126 molecular outflow. The values given in parentheses have been corrected for optical 
depth (assuming rco/(l ~ <?~ Tco ) = 3.5) and inclination (assuming i = 57.3°). 



Blue Red Total 



Maximum velocity (km/s) 


6 


(11.1) 


7 


(13) 






Size (AU/pc) 


20500/0.10 


(24400)/(0.12) 


19500/0.095 


(23200)/(0.11) 






Mass (xlCT 4 M ) 


1.6 


(5.5) 


8.1 


(28.2) 


9.6 


(33.7) 


Dynamical timescale (xlO 3 years) 


16.2 


(10.4) 


15.4 


(9.9) 






Momentum (xlO 3 M Q km/s) 


5.4 


(35) 


25.6 


(166) 


31.0 


(201) 


Momentum flux(xlO~ 8 M Q km/s y" 1 ) 


3.3 


(34) 


15.8 


(160) 


19.1 


(193) 


Kinetic Energy (xlO~ 3 M Q (km/s) 2 ) 


0.98 


(12) 


4.6 


(55) 


5.6 


(67) 


Mechanical luminosity (xKT 5 L Q ) 


1.0 


(19) 


4.9 


(92) 


5.9 


(110) 
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Fig. 13. CO J=3-2 velocity channel maps of the MMS126 out- 
flow. The greyscale is tailored to emphasize the low surface 
brightness outflow lobes. The contours are spaced by 1 K. 



7.1. The CO outflow 

We observed the area around MMS126 in the CO(3-2) line us- 
ing the JCMT 1 on June 17 & 30 and on July 30, 2005. We took 
two position switched on-the-fly maps centered on MMS126, 
each covering a 80"xl80" field at a position angle of 32°, sam- 
pling every 5" with a 5" spacing between rows. Smaller maps 

1 The James Clerk Maxwell Telescope is operated by The 
Joint Astronomy Centre on behalf of the Particle Physics and 
Astronomy Research Council of the United Kingdom, the Netherlands 
Organisation for Scientific Research, and the National Research 
Council of Canada. 



were added towards the end of the outflow lobes. The integra- 
tion time on each point was 4 sec. The weather was good, with a 
stable atmosphere and T230 ghz ~ 0.055 (June 17), ~ 0.08 (June 
30), and ~ 0.045 (July 30). In addition, spectra in frequency- 
switched mode were taken towards the position-switch OFF 
position. The spectra were baseline-subtracted, resampled to 
0.6 km/s resolution, the OFF-position spectrum added, and fi- 
nally gridded into one datacube within CLASS. 

Figure [O] shows the resulting channel maps. Along the 
outflow axis indicated by our near-infrared H2 observations 
there is clear evidence for a CO flow. Red-shifted emission is 
seen towards the south-west and some blue-shifted emission 
towards the north-east. Overall, the range over which higher- 
velocity gas is seen is fairly narrow. In the south-western, red- 
shifted lobe, it extends over a few velocity bins from about 5.7 
to 8.5 km/s. Reasonably bright emission in the north-eastern, 
blue-shifted lobe is only seen in its tip, otherwise it is hardly 
separated and only faintly seen in two channels around km/s; 
the contours in the velocity bins around 1 km/s indicate that 
this lobe contributes significantly at these velocities too, but is 
diluted by ambient emission. Outflow emission extends over 
about 2^5 (20000 AU/0. 1 pc) out to the edges of the map, with 
both lobes showing a somewhat brighter, compact CO blob 
there, which might well be a terminating working surface; H2 
emission is also seen out to a similar distance in the south- 
western lobe. 

We have derived estimates of basic flow properties as fol- 
lows: the measured brightness temperature was converted to 
molecular gas column densities assuming the CO to be opti- 
cally thin and in LTE at a temperature of 30 K; the CO abun- 
dance was assumed to be 10~ 4 of the H2 abundance. Masses 
were derived as a function of velocity for each 0.6 km/s wide 
channel. Kinetic energies and momenta were derived as E^ n = 
l/2£Af(v) • (v - v cen ) 2 and P = £Af(v) • (v - v ce „), with 
Veen = 3.4 km/s. A characteristic timescale was derived by di- 
viding the maximum distance at which high-velocity emission 
was seen by the maximum CO velocity observed. The south- 
western, red-shifted lobe is seen to extend to the edge of the 
mapped area, so a maximum distance of 90" was used; the 
north-eastern, blueshifted lobe seems to terminate within the 
mapped area. Finally, mechanical luminosities and momentum 
input rates were determined by dividing the total kinetic energy 
and momentum by the characteristic timescale. 

Table^gives a list of outflow properties, for the blue- and 
red-shifted lobes alone and for the entire outflow. In paren- 
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theses we also list the outflow parameters after applying cor- 
rections for optical d epth effects and inclination, following 
iBontemps et alJ ( fl99 6:) (assuming rco/(l -e~ Tco ) = 3.5 and an 
inclination i = 57.3°). The red- and blue-shifted emission does 
not seem to overlap, indicating that the flow axis is neither very 
close to the plane of the sky nor very nearly perpendicular to it, 
so i = 57.3° might in fact be a reasonable assumption. 

A comparison w ith the outflow parameters listed in 
IBontemps et alJ |l]?96) shows that the MMS 126 molecular out- 
flow is one of the less energetic flows, as could be anticipated 
by the poor contrast of the flow against the ambient cloud mate- 
rial even in the CO(3-2) line and its fairly narrow extent in ve- 
locity space. Its momentum flux rate is much more typical for a 
more evolved object than for a Class protostar. However, our 
map shows the flow to be fairly well collimated, which is more 
typical for the very youngest molecular outflows. Along with 
the deeply embedded, low circumstellar mass, low-luminosity 
nature of the driving source, the more likely explanation for 
the low energy/momentum input rate to the flow seems to be a 
very low protostellar mass (and age) rather than MMS 126 be- 
ing a relatively mature object at the end of the accretion and 
outflow phase. 

8. Summary & Conclusions 

We have performed a 1 .2 millimetre survey of the main body of 
the p Ophiuchi star-forming cloud. Our one square degree field 
enables us to verify and extend the results from previous sur- 
veys (MAN98 and J00). We detect and measure 143 cores, in- 
cluding 111 with no detected stellar objects. We identify a core, 
MMS 126, which appears to harbour a new low-mass Class 
protostar, driving a CO molecular outflow. 

By comparing core masses with their sizes (as measured by 
the flux and FWHMs), we find that the cores contained within 
an inner circle of 0.2° are compacter. This implies they are gen- 
erally denser and of higher pressure. However, on comparison 
to the surface pressure required for hydrostatic equilibrium, 
there appear to be vast pressure variations (by a factor 10 to 
100) within each location. 

The data indicate that the mean core mass increases with 
the core radius. This can be measured despite the expected 
lower bound to the core masses given by the observational sen- 
sitivity (which would tend to distort the relation towards the 
form log M oc log R e ) and despite the very wide range in possi- 
ble core masses at each radius. 

The clustering has also been investigated by deriving the 
two-point correlation function. We confirm previous results 
that the cores are highly clustered, similar to that of galaxy 
clusters (J00). However, from the two-point correlation func- 
tion, we find no evidence for a critical length to identify with 
a Jeans length. However, the frequency of nearest neighbours 
provides a more local measure and reveals an abundance of 
small cores with low separation (5000 AU) - of order of the 
mean core size. This suggests that larger cores may fragment 
on this scale. 

Moreover, the relevant velocity dispersion (in place of the 
sound speed) will also vary across the cloud, possibly produc- 
ing velocity variations positively correlated with separation. 



Therefore, although the Jeans length remains a relevant length 
scale in turbulent fragmentation and collapse theory, an ex- 
tended power-law correlation function could still be expected 
in a turbulent environment. 

The cores are not distributed isotropically. There is a 
strongly preferred direction, transverse to the direction of the 
Sco OB association. However, the orientations of the major 
axes of the individual cores are consistent with a random dis- 
tribution (|B] online version only). 

The mass function may be approximated by a broken power 
law. As previously noted by MAN98 for the inner regions, this 
is similar to that of the galactic field IMF for stars. The flatten- 
ing at lower masses occurs at a mass of ~0.1 to 0.3 M , about 
twice as high as the equivalent break mass found for stars. We 
suggest that this difference is consistent with the core mass 
loss expected during the star formation process although the 
absolute value of the break is sensitive to assumed mass con- 
version parameters. An upper break at ~0.5 - 1.0 M is also 
found. A smoothly changing function such as a log-normal 
might be better suited to describe the core mass function de- 
rived from our observations, similar to the mass functions de- 
rived from numerical models o f turbulence driven clouds by 
iBallesteros-Paredes et ail J2005I) . 

The core mass function does not vary with position. That is, 
there is no evidence for a general mass segregation although the 
few highest core masses are found in the inner crowded zone. 

How have the cores in L1688 developed? Turbulen t mo- 
tions are inherent to the entire clumps (Lor en et alJ 1 1 990) 
and (although subsonic) also to the cores feelk)che"etall 
2001). We may thus expect supersonic turbulence to drive 
high pressure and density variations on the scale of the cloud 
and clumps. Some density peaks within the chaotic flow 
may involve sufficient mass to become gravitationally bound; 
Elm egreen & S hadmehra J2003I) argue that this scenario is the 
most likely one of the three they consider to have produced 
the millimetre continuum sources in Ophiuchus. The turbu- 
lence then decays within these proto-cores although it is still 
sufficient to generate molecular line profil es and co ntribute to 
core support in the observed cores (Bell oche et aflbOOll) . As 
found in numerical simulations, the structure of these dynam- 
ical cores may closely resemble hydrostatic, gravitationally 
bound objects ( Ballest eros-Paredes et alJ 12003): they should 
also possess prominent infal l signatures as indeed are found in 
several starless cores teellocheet alfc oOl). Furthermore, it is 
very likely that the detected cores are biased towards the slow- 
est evolving long-lived objects. That is, the cores we actually 
observe are not typical of those which form but are simply a 
particular subset of cores which evolve very slowly, perhaps 
not even prone to collapse into stars. 

The present study is consistent with a supersonic turbulence 
interpretation through the wide range in pressures, the hierar- 
chical clustering and the randomness in core orientations, and 
the shape of the mass function. Turbulence, however, is gen- 
erally supposed to be a means of dynamical support and rapid 
dispersal rather than quiescent confinement of equilibrium core 
configurations. If the clouds are close to equilibrium, the lin- 
ear mass-radius relation would correspond to a fixed velocity 
dispersion. This indicates that a transition to coherence has al- 
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ready taken place and is a sign that the H2 velocity dispersion is 
transonic or subsonic in the cores ( Goodma n et alJl998l). In re- 
cent h ydrodynamic gravoturbulent simulations ( Kl essen et alJ 
120051) . the majority of cores are actually found to be coherent 
and 80% are subsonic or transonic. The cores are produced be- 
hind shock waves, at stagnation points in converging turbulent 
flows. However, starless cores are also found to be gravitation- 
ally unbound in these simulations. This should in the future be 
tested by measurement of the virial masses of a large sample of 
the starless cores. 

Further implications of the above results will be explored in 
a following work in which we relate the millimetre properties 
to those of the infrared and optical outflows from protostars 
and young stars. In combination with molecular spectroscopic 
studies of the cores, we will be able to constrain models for star 
formation in p Ophiuchus. 

Acknowledgements. TS thanks the Alexander von Humboldt 
Gesellschaft for support through a Feodor Lynen fellowship. MDS 
thanks DCAL, Northern Ireland. Thanks are due to Doug Johnstone 
for providing the J00 and J04 p-Oph SCUBA maps. This research has 
made use of the SIMBAD database, operated at CDS, Strasbourg, 
France. The Digitized Sky Survey was produced at the Space 
Telescope Science Institute under U.S. Government grant NAG 
W-2166. The images of these surveys are based on photographic 
data obtained using the Oschin Schmidt Telescope on Palomar 
Mountain and the UK Schmidt Telescope. The plates were processed 
into the present compressed digital form with the permission of 
these institutions. This publication makes use of data products from 
the Two Micron All Sky Survey, which is a joint project of the 
University of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National Science 
Foundation. 



References 

Allen, L. E., Myers, P. C, Di Francesco, J., et al. 2002, ApJ, 
566, 993 

Andre, P. & Montmerle, T. 1994, ApJ, 420, 837 
Andre, P., Ward-Thompson, D., & Barsony, M. 1993, ApJ, 406, 
122 

Andrews, S. M. & Williams, J. D. 2005, ApJ, 619, L175 
Ballesteros-Paredes, J., Gazol, A., Kim, J., & Klessen, R. S. 

2005, ApJ, in press (astroph/0509591) 
Ballesteros-Paredes, J., Klessen, R. S., & Vazquez-Semadeni, 

E. 2003, ApJ, 592, 188 
Barsony, M., Carlstrom, J. E., Burton, M. G., Russell, A. P. G., 

& Garden, R. 1989, ApJ, 346, L93 
Bate, M. R., Clarke, C. J., & McCaughrean, M. J. 1998, 

MNRAS, 297, 1163 
Belloche, A., Andre, P., & Motte, F. 2001, in ASP Conf. Ser. 

243: From Darkness to Light: Origin and Evolution of Young 

Stellar Clusters, 3 13-+ 
Bonnor, W. B. 1956, MNRAS, 116, 351 
Bontemps, S., Andre, P., Kaas, A. A., et al. 2001, A&A, 372, 

173 

Bontemps, S., Andre, P., Terebey, S., & Cabrit, S. 1996, A&A, 
311,858 



Brandner, W., Sheppard, S., Zinnecker, H., et al. 2000, A&A, 
364, L13 

Comeron, F, Rieke, G. H., Burrows, A., & Rieke, M. J. 1993, 
ApJ, 416, 185 

de Geus, E. J., de Zeeuw, P. T., & Lub, J. 1989, A&A, 216, 44 
de Zeeuw, P. T., Hoogerwerf, R., de Bruijne, J. H. J., Brown, 

A. G. A., & Blaauw, A. 1999, AJ, 117, 354 
Ebert, R. 1955, Zeitschrift fur Astrophysics, 37, 217 
Elias, J. H. 1978, ApJ, 224, 453 

Elmegreen, B. G. & Shadmehri, M. 2003, MNRAS, 338, 817 
Galli, D., Walmsley, M., & Goncalves, J. 2002, A&A, 394, 275 
Gomez, M., Hartmann, L., Kenyon, S. J., & Hewett, R. 1993, 
AJ, 105, 1927 

Goodman, A. A., Barranco, J. A., Wilner, D. J., & Heyer, M. H. 

1998, ApJ, 504, 223 
Grasdalen, G. L., Strom, K. M., & Strom, S. E. 1973, ApJ, 184, 

L53 

Johnstone, D., Di Francesco, J., & Kirk, H. 2004, ApJ, 611, 
L45 

Johnstone, D., Wilson, C. D., Moriarty-Schieven, G., et al. 

2000, ApJ, 545, 327 
Khanzadyan, T., Gredel, R., Smith, M. D., & Stanke, T. 2004, 

A&A, 426, 171 

Klessen, R. S., Ballesteros-Paredes, J., Vazquez-Semadeni, E., 

& Duran-Rojas, C. 2005, ApJ, 620, 786 
Kroupa, P. 2001, MNRAS, 322, 231 
Larson, R. B. 1995, MNRAS, 272, 213 
Loren, R. B. 1989a, ApJ, 338, 902 
Loren, R. B. 1989b, ApJ, 338, 925 

Loren, R. B., Wootten, A., & Wilking, B. A. 1990, ApJ, 365, 
269 

Luhman, K. L. & Rieke, G. H. 1999, ApJ, 525, 440 
Maddox, S. J., Efstathiou, G., Sutherland, W. J., & Loveday, J. 

1990, MNRAS, 242, 43P 
Mezger, P. G., Sievers, A., Zylka, R., et al. 1992, A&A, 265, 

743 

Motte, F, Andre, P., & Neri, R. 1998, A&A, 336, 150 
Nakajima, Y., Tachihara, K., Hanawa, T., & Nakano, M. 1998, 
ApJ, 497,721 

Niirnberger, D., Brandner, W., Yorke, H. W., & Zinnecker, H. 

1998, A&A, 330, 549 
Rebull, L. M., Wolff, S. C, & Strom, S. E. 2004, AJ, 127, 1029 
Reichertz, L. A., Weferling, B., Esch, W., & Kreysa, E. 2001, 

A&A, 379, 735 
Simon, M. 1997, ApJ, 482, L81 + 

Tachihara, K., Mizuno, A., & Fukui, Y. 2000, ApJ, 528, 817 
Vrba, F. J. 1977, AJ, 82, 198 

Vrba, F. J., Strom, K. M., Strom, S. E., & Grasdalen, G. L. 

1975, ApJ, 197, 77 
Weferling, B., Reichertz, L. A., Schmid-Burgk, J., & Kreysa, 

E. 2002, A&A, 383, 1088 
Wilking, B. A., Bontemps, S., Schuler, R. E., Greene, T. P., & 

Andre, P. 2001, ApJ, 551, 357 
Wilking, B. A., Lada, C. J., & Young, E. T. 1989, ApJ, 340, 

823 

Williams, J. P., de Geus, E. J., & Blitz, L. 1994, ApJ, 428, 693 
Wilson, C. D., Avery, L. W., Fich, M., et al. 1999, ApJ, 513, 
L139 



T. Stanke et al.: The signatures of protostars in p Ophiuchus II. Millimetre continuum observations, Online Material p 1 



Online Material 



T. Stanke et al.: The signatures of protostars in p Ophiuchus II. Millimetre continuum observations, Online Material p 2 

Appendix A: Notes on individual objects 




16 27 m 00 s 40 s 20 s 

Fig.A.2. The location, size, and orientation of the millimetre 
sources in and around the molecular clump A. 



Rho Oph B 




16 h 27 m 40 s 20 s 00 s 

Fig.A.3. The location, size, and orientation of the millimetre 
sources in and around the molecular clump B. 



In this section we present a list of sources identified from 
our survey (Tab. IA.lt . along with finding charts (Figs. IA.fllA.2l 
IA.3I and IA.4> and short notes on the properties of individual 
sources, such as the presence of optical or near-infrared (stel- 
lar) counterparts, or whether absorption features as seen on the 
DSS or 2MASS images can be associated. 



available as jpg: 0511093.1331fa04.jpg 

full postscript version available at 

http://www.ifa.hawaii.edU//users/stanke/preprints.html 



Fig. A.4. The location, size, and orientation of the millimetre 
sources in and around the molecular clumps C, E, and F. 

MMS001: ridge NE of VL A 1623: includes starless cores SMI 

and SM1N; no optical/IR stellar counte rpart. 

MMS002: VLA1623 Class jAndre et al.lll993h : no opt/ir 
counterpart. 

MMS003: pointlike source; optical/IR stellar counterpart 
YLW32. 

MMS004: SM2 starless core; no opt/ir counterpart 
MMS005: compact condensation in filament W of 
SM1/VLA1623 ridge; YSO LFAM1; no optical counter- 
part; bright NIR source, extended K-band nebula. 
MMS006: pointlike source; GSS39/Elia2-27; no opt. coutner- 
part; IR counterpart. 

MMS007: part of filament extending NE of SM1/VLA1623; 
no clear opt/IR stellar counterpart. 

MMS008: compact source: multiple T Tauri star SR24; the 
flux seems to be associated with the ( so far) unresolved south- 
ern single component, as claimed bv iNiirnb erger et aD {l998) 
and confirmed recently bv lAndrews & Williams J2005I) . rather 
than the northern (close binary) component. Our flux measure- 
ment i s in good agreement with the result of INiirnberger et all 
ill 9981) . 

MMS009: compact source; no opt. stellar counterpart; NIR 

star Elias 2-29 (also includes GY210). 

MMS010: compact source; optical/IR counterpart YLW34 

MMS011: compact source; YSO WLY 2-43; Class I; no opt. 

stellar counterpart; NIR star; high extinction. 

M MS012: compact, bu t clearly extended source; SMM 16267 

of IWilson et alJ ([1999) (suggesting it to be a prestellar core) 

no optical counterpart; high extinction; fuzzy K band emission 

south of clump 

MMS013: compact source; no opt. star; NIR star GSS26. 
MMS014: elongated feature SE of MMS005; no optical coun- 
terpart; some IR stars might be associated. Breaks up into chain 
of three subcores in MAN98. 

MMS015: bright compact core; no opt. stellar counterpart; 
NIR counterpart VSSG17; YSO; high opt. extinction. 
MMS016: bright elongated core; YSO VSSG18 appears to be 
associated with this feature. 

MMS017: part of filament extending NE of SM1/VLA1623; 
very faint K-band source at NE tip? 
MMS018: compact core; no opt./NIR stellar counterpart. 
MMS019: faint compact source; faint opt/bright NIR stellar 
counterpart YLW3 1 

MMS020: bright core; no clear opt./NIR stellar counterpart; 
high opt. extinction. 

MMS021: highly elongated; no opt./NIR stellar counterpart. 



Table A.l. List of 1.2 mm sources found in our survey. The sources are listed in the order they were identified by the wavelet decomposition plus clumpfind technique. F; nt 
gives the integrated source flux in mJy; F pea k (a), (b), (c), and (d) give the peak fluxes determined in the following ways: (a): the maximum pixel value inside an ellipse defined 
by the minor and major axes of a gauss-fit to the source; (b): as above, but using only the wavelet scales from where on the source is detected first (given in parenthesis); (c): 
the peak flux after re-binning according to the source minor axis (see Sect.[3]for details); (d): the peak flux of the composite Gaussian model for the source. The columns 
labelled maj and min give the major and minor axis of a single Gaussian fit to the source (in arcseconds). PA gives the posit ion angle of the source majo r axis in degrees east 
of north. The columns labelled as MAN98 and J00 give cross-identifications to the sources identified by Mot te et all Jl998i) and ljohnstone et all J2000I) . respectively. A long 
dash means that the source is outside the respective survey area, and "n.i." labels features which are visible in the respective survey, but not identified as a separate feature. 
Cross-identifications given in parentheses mean that a feature found here is a part of a feature noted by MAN98 or J00. Finally, a comment on association with known YSOs 
is given. Apparently starless features are marked as "si". 





RA 


DEC 


Pint 




F 

p l 


ak 


(mJy) 




FWHM (") 


PA 


MAN98 


J00 


Comment 




J2000 


J2000 


(mJy) 
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(c) 


(d) 


maj 


min 










MMS001 


16:26:27.2 


-24:23:46 


7775 


1999 


1999 


(1) 


1828 


2002 


65 


33 


-7 


SMI, SM1N 


16264-2423 


si (SMI, SM1N) 


MMS002 


16:26:26.2 


-24:24:28 


1205 


924 


924 


(1) 


864 


912 


29 


27 


85 


VLA1623 


16264-2424 


YSO: VLA1623 Class 


MMS003 


16:26:24.0 


-24:16:12 


242 


248 


248 


(1) 


227 


251 


26 


24 


-28 


— 


16264-2416 


YSO: YLW32 opt/IR TTS 


MMS004 


16:26:30.4 


-24:24:31 


10255 


973 


973 


(1) 


912 


987 


95 


57 


-63 


SM2 


n.i. 


si 


MMS005 


16:26:21.3 


-24:22:48 


1060 


357 


357 


(1) 


320 


351 


50 


35 


-29 


LFAM1 


16263-2422 


YSO: LFAM1; Class I? 


MMS006 


16:26:44.9 


-24:23:06 


291 


230 


230 


(1) 


217 


226 


28 


24 


6 


— 


16267-2423 


YSO: GSS39 IR star Class I? 


MMS007 


16:26:28.3 


-24:22:45 


3137 


490 


490 


(1) 


425 


482 


77 


51 


-41 


A-MM6/7 


n.i. 


si 


MMS008 


16:26:58.4 


-24:45:36 


228 


145 


145 


(1) 


133 


153 


30 


24 


59 


— 


16269-2445 


YSO: SR24 TTS binary 


MMS009 


16:27:09.1 


-24:37:25 


316 


154 


154 


(1) 


129 


158 


34 


30 


-1 


Elia 2-29 


16271-2437b 


YSO:Elia 2-29 


MMS010 


16:26:23.7 


-24:43:12 


153 


132 


132 


(1) 


119 


142 


26 


25 


-32 


— 


— 


YSO: YLW34 TTS 


MMS011 


16:27:27.3 


-24:40:55 


967 


193 


193 


(1) 


138 


193 


82 


58 


79 


IRS43 


16274-2440 


YSO: WLY 2-43 Class I 


MMS012 


16:26:43.8 


-24:17:22 


1109 


205 


205 


(1) 


170 


201 


73 


55 


43 


— 


16267-2417 


si 


MMS013 


16:26:10.3 


-24:20:52 


84 


70 


70 


(1) 


70 


78 


26 


19 


62 


GSS26 


16261-2420 


YSO: GSS26 Class I? 


MMS014 


16:26:22.5 


-24:23:48 


964 


265 


265 


(1) 


254 


261 


67 


32 


-19 


A-MM 1/2/3 


n.i. 


si 


MMS015 


16:27:29.6 


-24:27:37 


1336 


284 


284 


(1) 


246 


283 


56 


51 


-46 


B2-MM10 


16274-2427c 


YSO: VSSG17 (Class I?) 


MMS016 


16:27:27.8 


-24:27:13 


1058 


236 


236 


(1) 


220 


240 


75 


36 


63 


B2-MM8 


(16274-2427b) 


YSO: contains NIR YSO VSSG18 


MMS017 


16:26:30.8 


-24:22: 2 


1145 


187 


187 


(1) 


163 


202 


101 


43 


35 


n.i. 


n.i. 


si 


MMS018 


16:27:33.5 


-24:26:19 


2900 


350 


350 


u: 


307 


333 


78 


66 


-77 


B2-MM16/17 


16275-2426 


si 


MMS019 


16:26:18.7 


-24:28:18 


60 


76 


76 


(1) 


76 


76 


27 


21 


8 


— 


16263-2428 


YSO: YLW31 Class II? 


MMS020 


16:27:13.0 


-24:29:55 


1150 


152 


152 


(1) 


126 


155 


92 


66 


32 


B1-MM3 


(16272-2429) 


si 


MMS021 


16:27:33.3 


-24:27:32 


718 


130 


130 


(1) 


110 


125 


98 


44 


-9 


B2-MM15 


n.i. 


si 


MMS022 


16:26:58.6 


-24:34:06 


4181 


291 


291 


(1) 


242 


277 


112 


81 


38 


(C-S) 


(16269-2434) 


si 


MMS023 


16:27:16.1 


-24:31:13 


186 


83 


83 


(1) 


78 


78 


41 


31 


23 


B1-MM4 


(16272-2430) 


si 


MMS024 


16:26:23.2 


-24:24:43 


274 


127 


127 


(1) 


104 


101 


46 


35 


23 


LFAM3 


n.i. 


YSO: LFAM3 Class I ? 


MMS025 


16:26:43.3 


-24:34:50 


415 


117 


117 


(1) 


74 


117 


109 


68 


74 


WL12 


16267-2434 


YSO: WL12 


MMS026 


16:26:23.2 


-24:22:04 


601 


170 


170 


(1) 


161 


176 


72 


28 


25 


A-MM4 


16264-2421 


si 


MMS027 


16:27:15.4 


-24:30:34 


1651 


218 


218 


1) 


190 


207 


73 


72 


-39 


(Bl) 


(16272-2430) 


si 


MMS028 


16:27:11.6 


-24:29:15 


874 


176 


176 


(1) 


151 


161 


67 


50 


2 


B1-MM2 


(16272-2429) 


si 


MMS029 


16:27:28.9 


-24:26:39 


610 


225 


225 


(1) 


208 


238 


49 


30 


76 


B2-MM9/12 


(16274-2427b) 


si 


MMS030 


16:26:08.0 


-24:20:20 


2614 


146 


146 


(1) 


111 


153 


142 


92 


10 


Oph-A3 


16261-2419 


si 


MMS031 


16:27:24.2 


-24:27:42 


1318 


194 


194 


(1) 


174 


189 


88 


49 


39 


B2-MM3/4/5 


(16274-2427a) 


si 


MMS032 


16:26:40.4 


-24:27:13 


114 


83 


83 


(1) 


72 


85 


33 


28 


-58 


CRBR42 


16266-2427 


YSO: CRBR42 Class I? 


MMS033 


16:27:20.1 


-24:27:09 


624 


138 


138 


(1) 


124 


123 


76 


38 


45 


B2-MM2 


16273-2427 


si 


MMS034 


16:26:17.0 


-24:23:40 


126 


70 


70 


(1) 


67 


69 


26 


24 


-7 


CRBR12 


n.i. 


YSO: CRBR 12 Class I? 


MMS035 


16:27:23.8 


-24:40:54 


921 


167 


167 


(1) 


136 


161 


65 


64 


36 


CRBR85 + F-MM2 


n.i. 


YSO: Class I (?) 



Table A.2. List of 1.2 mm sources found in our survey. 
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Table A.3. List of 1 .2 mm sources found in our survey. 





RA 


DEC 


Pint 




Fpeak 


(mJy) 




FWHM (") 


PA 


MAN98 


J00 


Comment 




J2000 


J2000 


(mJy) 
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(b) 


(c) 


(d) 


maj 


min 










MMS075 


16:29:01.4 


-24:21:11 


37 


38 


25 (2) 


38 


45 


42 


18 


14 






si 


MMS076 


16:25:56.3 


-24:20:48 


102 


42 


30 (2) 


34 


38 


33 


27 


37 




16259-2420 


YSO: TTS HBC259 Class II? 


MMS077 


16:26:49.0 


-24:26:55 
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48 


38(2) 


43 
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77 


41 
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n.i. 


si 
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16:27:39.9 


-24:29:26 
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24(2) 


35 


37 
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26 
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n.i. 


n.i. 


si 
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-24:18:17 


355 


73 


56 (2) 


60 


62 


93 


43 
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n.i. 


si 
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16:27:30.5 


-24:41:56 


779 


77 


67 (2) 


60 


80 
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80 
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n.i. 


n.i. 


si 
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16:27:48.3 
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51 (2) 


45 


55 
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87 
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si 
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47 


34 (2) 


43 


40 


167 


33 


69 


-E-MMl 


n.i. 


si 


MMS083 


16:27:38.2 


-24:27:38 


180 


59 


43 (2) 


50 


58 


93 


34 


90 


n.i. 


n.i. 


si 


MMS084 


16:27:37.8 


-24:46:05 


894 


63 


44(2) 


40 


46 


200 


85 


-74 


— 


n.i. 


si 


MMS085 


16:26:30.4 


-24:18:39 


406 


69 


56 (2) 


59 


55 


75 


52 


87 


— 


n.i. 


si 


MMS086 


16:27:19.1 


-24:31:26 


118 


62 


46 (2) 


56 


60 


55 


29 


-7 


(Bl) 


(16272-2430) 


si 


MMS087 


16:26:47.1 


-24:29:28 


1246 


98 


74 (3) 


87 


85 


125 


71 


-39 


Oph-AC2 


16267-2429 


si 


MMS088 


16:25:58.6 


-24:18:47 


1786 


88 


64(3) 


69 


65 


259 


90 


-57 


n.i. 


n.i. 


si 


MMS089 


16:27:25.6 


-24:33:32 


640 


56 


38(3) 


44 


46 


90 


73 


-77 


— 


n.i. 


si 


MMS090 


16:27:01.8 


-24:26:25 


720 


66 


41 (3) 


46 


56 


136 


72 


-89 


n.i. 


n.i. 


si 


MMS091 


16:25:58.8 


-24:31:29 


1392 


81 


45 (3) 


41 


41 


178 


146 


-29 


— 


n.i. 


si 


MMS092 


16:29:19.5 


-24:41:33 


504 


62 


29 (3) 


33 


44 


85 


67 


-46 


— 


— 


si 


MMS093 


16:27:43.6 


-24:13:32 


1759 


74 


46 (3) 


44 


46 


176 


111 


-53 


— 


— 


si 


MMS094 


16:26:59.0 


-24:22:59 


1584 


64 


50 (3) 


53 


52 


222 


94 


67 


Oph-AB 


n.i. 


si 


MMS095 


16:27:33.3 


-24:15:40 


992 


56 


37 (3) 


38 


46 


177 


115 


-27 


— 


— 


si 


MMS096 


16:26:55.4 


-24:19:14 


545 


51 


34 (3) 


38 


37 


146 


90 


42 


Oph-G 


16269-2419 


si 


MMS097 


16:28:14.2 


-25:02:58 


589 


70 


29 (3) 


41 


31 


138 


69 


-19 


— 


— 


si 


MMS098 


16:25:47.3 


-24:29:54 


573 


69 


26 (3) 


34 


25 


201 


85 


66 


— 


n.i. 


si 


MMS099 


16:28:07.4 


-24:35:20 


2658 


64 


46 (3) 


45 


55 


283 


125 


-64 


— 


— 


si 


MMS100 


16:27:56.4 


-24:42:48 


637 


47 


28 (3) 


36 


34 


212 


64 


82 


— 


— 


si 


MMS101 


16:27:12.7 


-24:21:23 


648 


49 


32 (3) 


34 


34 


127 


96 


-88 


— 


— 


si 


MMS102 


16:29:07.7 


-24:35:03 


324 


50 


21 (3) 


41 


25 


135 


48 


-88 


— 


— 


si 


MMS103 


16:27:55.2 


-24:35:07 


323 


50 


29 (3) 


44 


38 


86 


51 


-75 


— 


— 


si 


MMS104 


16:26:27.9 


-24:33:36 


1166 


65 


41 (3) 


42 


46 


141 


119 


-54 


— 


n.i. 


si 


MMS105 


16:25:44.6 


-24:16:56 


1068 


55 


32 (3) 


32 


31 


191 


131 


49 


— 


n.i. 


si 


MMS106 


16:26:48.9 


-24:23:33 


169 


39 


25 (3) 


36 


38 


64 


50 


49 


n.i. 


n.i. 


si 


MMS107 


16:27:40.3 


-24:35:50 


1488 


57 


39 (3) 


44 


40 


277 


104 


-43 


— 


n.i. 


si 


MMS108 


16:27:20.6 


-24:24:03 


521 


50 


32 (3) 


42 


42 


123 


64 


-19 


(B3) 


n.i. 


si 


MMS109 


16:29:11.2 


-24:40:58 


590 


47 


28 (3) 


33 


28 


157 


98 


-75 


— 


— 


si 


MMS110 


16:25:49.4 


-24:09:34 


255 


40 


20(3) 


26 


26 


176 


64 


66 






si 


MMS111 


16:27:26.1 


-24:36:21 


180 


34 


19(3) 


29 


26 


111 


53 


68 




n.i. 


si 


MMS112 


16:27:10.7 


-24:26:15 


273 


49 


28 (3) 


40 


38 


95 


50 


28 


n.i. 


n.i. 


si 


MMS113 


16:26:34.3 


-24:35:27 


341 


51 


24(3) 


32 


26 


115 


68 


-58 




n.i. 


si 


MMS114 


16:27:28.0 


-24:20:32 


223 


40 


21 (3) 


29 


26 


88 


58 


59 






si 


MMS115 


16:26:16.2 


-24:31:21 


757 


58 


30 (3) 


31 


31 


158 


109 


-76 




n.i. 


si 



Table A.4. List of 1.2 mm sources found in our survey. 





RA 


DEC 


Fint 




Fpeak (mJy) 




FWHM (") 


PA 


MAN98 


J00 


Comment 




J2000 


J2000 


(mJy) 


(a) 


(b) 


(c) 


(d) 


maj 


min 










MMS116 


16:29:01.7 


-25:06:30 


1214 


76 


37 (3) 


38 


39 


181 


132 


-67 


— 


— 


si 


MMS117 


16:27:39.1 


-24:41:15 


315 


49 


31(3) 


43 


43 


71 


54 
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Fig. A.l. The location, size, and orientation of the 143 millimetre sources identified from our survey. 



MMS022: bright large core; no opt./NIR stellar counterpart; 
high extinction. 

MMS023: compact source; no clear optical counterpart; possi- 
bly IR star associated; high opt. extinction. 
MMS024: LFAM3 faint opt., bright nebulous NIR source; vis- 
ible on J00 map, but not noted as separate feature 
MMS025: compact source (possibly surrounded by extended 
Halo); no opt. stellar counterpart; IR stellar counterpart WL12; 
high opt. extinction. 

MMS026: elongated core NE of MMS005; no optical counter- 
part; faint K-band source? 

MMS027: bright core; no clear opt./NIR stellar counterpart; 
possibly IR star to the SE associated; high opt. extinction. 
MMS028: compact core; no clear opt./NIR stellar counterpart; 
high opt. extinction. 

MMS029: elongated core. No opt./NIR stellar counterpart. 
MMS030: extended core at SE tip of ridge (MMS088); no re- 
lated optical feature. 

MMS031: elongated core; no opt./NIR stellar counterpart. 
MMS032: faint compact source; no optical counterpart; faint 
IR star CRBR42; ISO-Oph 54; [GY92] 9 1 
MMS033: elongated core; no opt./NIR stellar counterpart. 
MMS034: faint compact source; no optical counterpart; IR 
stellar counterpart CRBR12. 

MMS035: compact source; no opt. stellar counterpart; NIR 
stellar counterpart, YSO CRBR85 Class I?; splits into two 
cores in MAN98 (YSO CRBR85 and F-MM2 starless core). 
MMS036: no stellar counterpart for main core; MAN98 E- 
MM2b might have faint K-band counterpart. 
MMS037: elongated core; no opt./NIR stellar counterpart. 
MMS038: bright compact core; no opt./NIR stellar counter- 
part; high extinction. 

MMS039: bright core adjacent to MMS022; no opt./NIR stel- 
lar counterpart; high extinction. 

MMS040: elongated core; no opt./NIR stellar counterpart; 
high extinction. Contains MAN98 starless core F-MM1 and, 
at its north-western tip, the YSO CRBR72 (Class I). 
MMS041: bright core; no opt./NIR stellar counterpart; near tip 
of cometary shaped high-optical-extinction cloud. The edge of 
the cloud are outlined by faint extensions to the east and west 
of this core. 

MMS042: compact source; YSO LFAM26; no opt. stellar 
counterpart; faint K-Band source; optical extinction ridge. 
MMS043: compact faint source, YSO SR21 
MMS044: large core; compact source at center? no opt./NIR 
stellar counterpart; high extinction. 

MMS045: faint compact core; no opt. stars; 3 NIR YSO stellar 
sources: IRS44/IRS46; CRBR88; 3 subsources det. by MAN98 



MMS046: faint compact source; no clear opt/NIR counterpart. 
MMS047: bright core; no opt./NIR stellar counterpart; at tip of 
cometary dark cloud. 

MMS048: small core; no optical/NIR counterpart; high ex- 
tinction. Part of small clump (together with MMS079 and 
MMS085). 

MMS049: small faint core; no opt./NIR stellar counterpart; 
high opt. extinction. 

MMS050: core in SMI ridge; no clear opt/NIR stellar counter- 
part. 

MMS051: compact core; no opt./NIR stellar counterpart; high 
extinction. 

MMS052: rather bright elongated core; no opt./NIR stellar 
counterpart; tail of cometary dark cloud. 
MMS053: small faint core; high optical extinction; no stellar 
opt./NIR counterpart. 

MMS054: faint core; high optical extinction; no clear opt/IR 
stellar counterpart. 

MMS055: faint core; no clear opt./NIR stars associated; high 
extinction. 

MMS056: diffuse core; no opt./NIR stellar counterpart; high 
opt. extinction. 

MMS057: small core; no clear optical/NIR counterpart. Part of 
small core (together with MMS063 and MMS069). 
MMS058: faint compact source; YSO WSB52/LFAM 
p8/[GY92] 314; opt./NIR star. 

MMS059: no opt./NIR stellar counterpart; high opt. extinction. 
MMS060: faint compact source; opt/NIR stellar counterpart 
WSB 60 = BKLT J162816-243657. 

MMS061: elongated core; superposition of several NIR YSOs: 
WL3/WL4/WL5/IRS37; no optical stellar counterpart; high ex- 
tinction. Class II sources. MAN98 classify this feature as com- 
posite pre-collapse core, but it might actually just be a super- 
position of circumstellar dust emission from the 4 NIR YSOs 
visible in this area. 

MMS062: compact source; associated with edge-on disk YSO 
jBrandner et alJ2000l) . 

MMS063: small core; no clear optical/NIR counterpart. Part of 
small clump (together with MMS057 and MMS069). 
MMS064: faint compact source; no opt. star; NIR stellar coun- 
terpart; YSO: IRS51 = WLY 2-51 = YLW 45. 
MMS065: faint compact source; seems to be associated with 
IR RN; YSO? hardly visible on MAN98 map; visible on J00 
map, but not noted as separate feature. 

MMS066: faint elongated core; no opt./NIR stellar counter- 
part; high opt. extinction. 
MMS067: no opt./NIR stellar counterpart. 
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MMS068: extended core; no clear opt./IR stellar conterpart. 
There is emission visible in the MAN98 map, but there is some 
morphological difference from our and the J00 maps. 
MMS069: small core; no clear optical/NIR counterpart. Part of 
small clump (together with MMS057 and MMS063). 
MMS070: faint core; core + envelope structure? No opt. stellar 
counterpart; NIR star WLY 2-54. 

MMS071: elongated core; no opt./NIR stellar counterpart; 
high extinction. Might split into a chain of faint sources. 
Includes NIR YSO WL20 at its south-eastern tip. 
MMS072: faint core at eastern end of clump B2; no opt./NIR 
stellar counterpart. 

MMS073: elongated core; no direct stellar counterpart, but 
[GY92] 301 is within core area. 

MMS074: filament extending north of SMI core; GSS31 on 
base of filament; possibly very faint nebulosity at tip. 
MMS075: faint compact source nearby MMS041; no opt./NIR 
stellar counterpart. 

MMS076: faint compact source; bright opt/NIR counterpart; 
V* V2058 Oph; HBC 259; SR4 

MMS077: small faint core; located on high-extinction fila- 
ment; no IR counterpart; not clearly visible in MAN98 map, 
part of their clump Oph-AC2 

MMS078: faint small core, no opt./NIR counterpart. 
MMS079: small core; no optical/NIR counterpart; high ex- 
tinction. Part of small clump (together with MMS048 and 
MMS085). 

MMS080: faint large core; no opt./NIR stellar counterpart. 
MMS081: large, elongated core; 3 NIR (YSO?) stars on our 
around the core: GY344, GY370, GY352; high opt. extinction. 
MMS082: faint filament; no opt/NIR stellar counterpart; high 
extinction. Corresponds roughly to MAN98 E-MM1, but dif- 
fers in morphology from our and the J00 maps. 
MMS083: faint core at eastern end of clump B2; no opt./NIR 
stellar counterpart. 

MMS084: large core; no clear opt./NIR stellar counterpart; 
close to edge of high-extinction cloud; close to edge of J00 
map; there is some emission, but it is not marked as separate 
feature. 

MMS085: small core; no optical/NIR counterpart; high ex- 
tinction. Part of small clump (together with MMS048 and 
MMS079). 

MMS086: compact source; no clear opt./NIR counterpart; high 
opt. extinction. 

MMS087: large core; no opt./NIR stellar counterpart; high ex- 
tinction patch. 

MMS088: SE-NW elongated ridge bounded by MMS030 and 
MMS105; apparently illuminated from the SW by HD147889 
(SRI) 

MMS089: faint extended core; no clear opt./NIR stellar coun- 
terpart; some bright NIR stars around; extinction patch. 
MMS090: faint core; on optical extinction filament; no clear 
opt/IR stellar counterpart. 

MMS091: faint large core; no stellar counterpart; optical ex- 
tinction patch. 

MMS092: faint core, no opt./NIR stellar counterpart; 
MMS092, MMS109, and MMS122 form a filament associated 
with high extinction. 



MMS093: low surface brightness core; on high optical extinc- 
tion lane (along with MMS095, MMS125, and MMS137); no 
clear opt/NIR stellar counterpart. 

MMS094: part of larger ridge; no opt./IR stellar counterpart; 
on opt. extinction lane. 

MMS095: low surface brightness core; on high optical extinc- 
tion lane (along with MMS093, MMS125, and MMS137); no 
clear opt/NIR stellar counterpart. 

MMS096: small low surface brightness core; no clearly asso- 
ciated optical/NIR counterpart. 

MMS097: extended low surface brightness feature; no clear 
feature seen on DSS/2MASS. 

MMS098: faint core; no stellar counterpart; optical extinction 
patch. 

MMS099: ridgelike, elongated core; no opt./NIR stellar coun- 
terpart; high-extinction lane. 

MMS100: large elongated core; no opt./NIR stellar counter- 
part; high extinction. 

MMS101: part of larger ridge; on high optical extinction lane; 
no clear opt/NIR stellar counterpart; on edge of J00 map. 
MMS102: very faint large core; no opt./NIR stellar counter- 
part; extinction patch. 

MMS103: faint core; no opt./NIR stellar counterpart; high ex- 
tinction. 

MMS104: faint large core; no opt/NIR stellar counterpart; 
high-extinction. 

MMS105: extended low surface brightness feature at NW end 
of ridge (MMS088); associated with diffuse K-band RN at SW 
edge; maybe some optical extiction patches associated. 
MMS106: core on larger ridge; no opt./IR stellar counterpart; 
on opt. extinction lane. 

MMS107: faint core; no clear opt./NIR stellar counterpart, but 
there are some IR stars in the core area which might belong to 
this feature; optical extinction filament. 

MMS108/MMS118/MMS143: faint cores; no opt./NIR stellar 
counterpart; high optical extinction, part of MAN98 B3. 
MMS109: faint core, no opt./NIR stellar counterpart; 
MMS092, MMS109, and MMS122 form a filament associated 
with high extinction. 

MMS110: small faint core; small optical extinction whisps; IR 
star near peak position. 

MMS111: faint core; no opt./NIR stellar counterpart; high ex- 
tinction. 

MMS112: faint core; no clear opt./NIR stellar counterpart; 
high opt. extinction lane? 

MMS113: faint large core; possibly faint compact source also; 
no opt./NIR stellar counterpart; high-extinction lane. 
MMS114: faint core; located on high optical extinction fila- 
ment; no clear opt./NIR stellar counterpart. 
MMS115: faint, large core; no opt./NIR stellar counterpart; 
high extinction lane. 

MMS116/120: extended low surface brightness features; no 
clear features seen on DSS/2MASS. 

MMS117: faint large core; no opt./NIR stellar counterpart; 
high extinction. 
MMS118: see MMS108 

MMS119: very faint large core; no opt./NIR stellar counter- 
part; extinction patch. 
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MMS120: see MMS116 

MMS121: faint core; no stellar counterpart, roughly corre- 
sponds to optical extinction patch. 

MMS122: faint core, no opt./NIR stellar counterpart; 
MMS092, MMS109, and MMS122 form a filament associated 
with high extinction. 

MMS123: faint core; no optical counterpart; NIR star within 
core area; high extinction. Faint on J00 map. 
MMS124: faint core, no opt./NIR stellar counterpart; high ex- 
tinction. 

MMS125: faint core; on high optical extinction lane (along 
with MMS093, MMS095, and MMS137); no clear opt/NIR 
stellar counterpart. 

MMS126: compact source; no opt./NIR stellar counterpart; 
IRAS 16253-2429? Class 0? Main source might have a com- 
panion to the north-east. 

MMS127: faint core; the YSO YLW37 (BKLT J162646- 
241203) is located at its north-eastern periphery. 
MMS128: very faint large core; no opt./NIR stellar counter- 
part; extinction patch. 

MMS129: very faint large core; no opt./NIR stellar counter- 
part; extinction patch. 

MMS130: large, faint, elongated core; no opt./NIR star associ- 
ated; extinction lane. 

MMS131/132/133/139: extended low surface brightness fea- 
tures; no clear features seen on DSS/2MASS. 
MMS134: low surface brightness; no opt. stellar counterpart, 
diffuse nebulosity; NIR star GSS15/BKLT J162535-243400 at 
north-eastern periphery. Faint stuff on J00 map. 
MMS135: low surface brightness; small optical extinction 
patch. 

MMS136: low surface brightness; 2 optical/N IR stars, V 2059 
Oph and ISO-Oph 180 jBontemps et alJ 2001) within core area. 
On edge of MAN98 and J00 maps. 

MMS137: low surface brightness; on high optical extinction 
lane (along with MMS093, MMS095, and MMS125); no clear 
opt/NIR stellar counterpart. 

MMS138: low surface brightness; some opt./NIR stars visible 
at core periphery. 
MMS139: see MMS131 

MMS140: faint compact source plus more extended core; no 
opt. counterpart; NIR star WL17; high extinction. 
MMS141: faint compact source; stellar NIR counterpart; opti- 
cally very faint; Class I/II? 
MMS142: faint extended core. 
MMS143: see MMS108 

Appendix B: Orientation 

On the largest scale, the p Ophiuchi cloud possesses southeast- 
northwest ridges (general orientation of ~ 145°) along with fil- 
amentary tails stretching in the NN E direction (~ 70°) over four 
degrees on the sky JLorenll989al) . 

A preferred cloud direction on sub-degree scales is also 
apparent, as shown on our 1.2 mm image (Fig.Q. Here, we 
quantify this by determining the orientation between each pair 
of cores. The resulting number distribution is displayed in 
Fig.lB~T1 



o 

Q_ 



CD 

E 

3 



1000 



800 



600 



400 



200 




50 100 150 

Pair Orientation (degrees) 

Fig. B.l. The number distribution of core-core orientations 
determined for each pair of cores (a total of 10,153 pairs; 
solid line), each pair separated by less than 30,000 AU (508 
pairs; dotted line) and those separated by more than 30,000 AU 
(dashed line). The latter two distributions have been normalised 
so that the total number is constant (i.e. displaying equal areas 
under the lines). 

Not surprisingly, the result is that the preferred core-pair 
direction is not only maintained on the sub-degree scale (solid 
line) but remains just as prominent below the separation of 
30,000 AU (0.06°). Thus the preferred direction of ~ 130 - 150° 
corresponds to the ridge direction - the expected direction of 
compression from shocks transmitted from the Sco OB2 asso- 
ciation JVrball977tlLorenll989bl) . 

Is the preferred direction still present on even smaller 
scales? To answer this, we plot the number distribution of the 
orientations of the major axis of each core in Fig. IB.2l We find 
no obvious trend in the data although the statistics are quite 
small. This result applies to the entire sample as well as the 
starless cores which implies that any preferred direction is lost 
at any early stage in core development. 

To further check for the presence of local (or global) align- 
ment of core orientations, we have analysed the absolute dif- 
ference in position angles (APA) that every core includes with 
its neighbours. Fig. IB. 31 (upper panel) shows the average po- 
sition angle difference of all cores with respect to their n'th 
nearest neighbour. The large scatter in APA around 45° indi- 
cates that the large majority of cores is randomly oriented with 
respect to their neighbours. At most a tiny tendency towards 
lower average APA is seen for the few very next neighbours; 
the still large scatter indicates that this can only be due to very 
few cases showing some degree of alignment, while the great 
majority is randomly oriented. 

Fig. lB.3l (lower panel) shows the average APA as a function 
of the projected distance from every core. Again, at most a very 
small trend towards smaller average APA at the smallest sepa- 
rations can be seen. The large scatter again indicates that this 
trend is caused by only a few core pairs, with most core pairs 
including random APA. 
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Fig. B.2. The number distribution of orientations of the major 
axes of the 111 starless cores (upper panel) and the entire 143 
cores (lower panel). 

Thus it appears that there is virtually no tendency for align- 
ment of core orientations, supporting the idea of core formation 
as a result of unordered, turbulent motions. In particular, sup- 
posing that the orientation of the core major axis is perpendicu- 
lar to the direction of preferred contraction, this contradicts the 
idea of ordered magnetic fields as producing agent of aligned 
protostellar disk/jet systems. 

Two (counteracting) effects should be kept in mind in the 
interpretation of this result: residual scanning effects might lo- 
cally push measured core orientations towards a similar posi- 
tion angle. On the other hand, we only can analyse the projected 
core distribution; a core seen to be the nearest neighbour in this 
2-D projection might in fact be at much larger distance in 3- 
D, which will tend to dilute any local alignment, thus in fact 
alignment might be more prevalent than seen here. 
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Fig. B.3. The relative orientation (difference in position angles) 
of core major axes with respect to their neighbours: average 
relative orientation with respect to the nth neighbour (upper 
panel), and average relative orientation of cores within a certain 
(projected) distance to the core. The smallest cores (with major 
axes determined to be smaller than 28") have been omitted. 
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